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PREFACE 



This publication contains the proceedings of a workshop conducted at the NASA 
Langley Research Center on the technology for controlling large space structures cur- 
rently being conceived for space applications by the year 2000. The purpose of the 
program was to address the fundamental technolo gy defi ciencies that were identified 
in several studies on large space ^t^~(LSSJ conducted by NASA irT the l^s"t Several 
years. Two_such studies haye_ pointed out the need for distributed control technology 
and technology for adaptive surface control. During the course of the workshop, addi- 
tional technological deficiencies were uncovered, the most notable being the lack of 
experience in managing large numbers of control system components and maintaining 
system level functions in the presence of failed actuator or sensor components. At 
present, these issues are of primary concern because there is little or no experience 
with them m the space program but the missions studies are critically dependent upon 
these technology items for their mission economy or success. 

A basic research program has been assembled at NASA Langley Research Center to 
address the technology deficiencies discussed previously. It consists of an in-house 
effort, university grants, and industry contracts. The staffs of the respective par- 
ticipants were assembled at the workshop to review the current state of research in 
the control technology for large structural systems and to plan the efforts that would 
be pursued by their respective organizations. This document contains the more impor- 
tant slides that were used by each participant with a word description where required 
for clarity. It is our intention to review the progress of the activity and have 
another workshop to plan the program for the following year at the close of the 1 981 
fiscal year (about October 1981). 

Use of trade names or names of manufacturers in this report does not consti- 
tute an official endorsement of such products or manufacturers, either expressed or 
implied, by the National Aeronautics and Space Administration. 

Raymond C. Montgomery, Cochairman u^ii 
Garnett C. Horner, Cochairman 
E. Burton Lightner, Coordinator 
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LANGLEY RESEARCH CENTER FLEXIBLE BEAM EXPERIMENT 



Raymond C. Montgomery 
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October 30-31, 1980 



LaRC FLEXIBLE BEAM EXPERIMENT 

To verify some of the optimization results and other control algorithms, 
a flexible beam experiment has been initiated at LaRC. In figure 1, the 
flexible beam experiment consists of a 3.66 m (12 ft) long aluminum beam with 
a 4.76 mm (3.16 in.) by 15 cm (6 in.) cross section. The beam is suspended by 
two small flexible cables so that free-free end conditions are approximated . 
Located in front of the beam are four electromagnetic shakers (actuators) 
which can be repositioned along the beam by sliding them along the platform 
which supports them. The console on the left contains the power amplifiers 
for the shakers. 





/ 



Figure 1 



LaRC FLEXIBLE BEAM EXPERIMENT 
(Continued) 

Figure 2 shows another picture of the experimental setup. On one side 
of the beam the four shakers are located and on the other side of the beam 
there are nine noncontacting displacement probes. With the experiment being 
tied in with the CDC Cyber 175 computer, real-time calculations may be made. 
For example, the output of the displacement probes can be made available to the 
computer. Using state estimation, the velocity at the shaker locations can 
be approximated. Knowing the damping rate or gain from the optimization pro- 
gram and the velocity, the desired force output of the shakers can be calcu- 
lated. 




Figure 2 
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G. C. Horner 
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OPTIMUM DAMPER LOCATIONS FOR A FREE- FREE BEAM 



The objectives of this research are to identify optimum locations for 
sensors and actuators on large space structures. If it is assumed that large 
platforms and antennae will have many potential actuator/sensor locations, we 
may logically ask "Where should actuators and sensors be placed?" Not only 
should the optimum placement be determined, but also the dynamic characteristics 
of actuators may also be necessary. 



OBJECTIVES 

o DEVELOP ALGORITHMS TO OPTIMALLY LOCATE AND DESIGN DAMPERS FOR LARGE 
SPACE STRUCTURES 

o DETERMINE REQUIREMENTS FOR DISTRIBUTED SENSING AND ACTUATION (AS 
OPPOSED TO COLOCATED SENSOR AND ACTUATOR) IN CONTROL OF STRUCTURAL 
SYSTEMS 

APPROACH 

o USE MATHEMATICAL PROGRAMMING TO SOLVE FOR OPTIMUM DAMPING RATE AND 
LOCATION, 

o CONSIDER ACTUATOR DYNAMICS TO SOLVE FOR OPTIMUM ACTUATOR MASS. 

Figure 1 



DAMPING CHARACTERISTICS OF A FREE- FREE BEAM 



To get an understanding of the behavior of large space structures, we 
first look at the damping characteristics of a uniform beam. A dash pot is 
located at one end of a free-free beam. This is an ideal dash pot which is 
characterized by a damping rate, C, and no other dynamic characteristics. In 
figure 2 it is seen that for small values of C (<.005), the damping ratio, % , 
and damping rate are linearily related. This is denoted as perturbation 
theory. As the damping rate is increased, the damping ratio reaches a peak 
value and then decreases. The peak value of the damping ratio is about 0.2 
for the first flexible mode. Suppose a design problem were stated which re- 
quired that the first mode have a damping ratio greater than 0.2. This 
requirement may be a result of mission performance specifications. To achieve 
more than the 0.2 damping ratio in the first mode, one or more dash pots are 
required. Since the design problem being addressed here is one in which the 
damping ratio is prescribed for each mode to be damped, the damping rate of 
the dash pots is determined. 
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DAMPING CHARACTERISTICS OF A CLAMPED- FREE BEAM 

The results are essentially the same as for the free-free beam in 
figure 2. 
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Figure 3 



NONLINEAR OPTIMIZATION PROBLEM 



A. design problem is posed which states that given the prescribed modal 
damping ratio for N modes, what are the optimum damping locations and sizes? 
The design problem is now cast as "a nonlinear optimization problem. Since it 
is not known where the dash pots should be located on a structure, the initial 
step is to put a dash pot at every location of the beam. The objective function 
is to minimize the total dissipative effort. The constraints are that the 
actual computed modal damping ratios must be greater than or equal to the 
prescribed value. Another constraint is that the damping rate must be 
positive. This guarantees stability. 



FOR PRESCRIBED MODAL DAMPING RATIO IN N MODES, WHAT ARE THE BEST DAMPING SIZES 
AND LOCATIONS? 
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Figure 4 



OPTIMUM DAMPING LOCATIONS AND SIZES FOR A FREE- FREE BEAM 



Some results are presented in figure 5 for a free-free beam. The design 
problem consisted of prescribing a modal damping ratio of 0.5 in N modes. The 
results are shown for N = 1, 2, 3, 4. The results are also split between 
symmetric solutions and nonsymmetric solutions. The symmetric solutions are 
obtained by minimizing the total dissipation while imposing symmetry in the solu- 
tion. The horizontal lines represent the length of the beam. The vertical 
lines are proportional to the magnitude of the damping rate at the location 
shown on the beam axis. The nonsymmetric solution is obtained by removing the 
symmetry requirement and the smallest damper location. Thus, nonsymmetric solu- 
tions will have no more than one fewer dampers than the symmetric case. In 
some cases the objective function for the nonsymmetric solution is less than 
that for the symmetric case. 
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OPTIMUM DAMPER LOCATIONS AND SIZE FOR A CLAMPED- FREE BEAM 



The results shown in figure 6 are similar to those in figure 5. 
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Figure 6 
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PUTURE RESEARCH 



The future research thrusts will involve the addition of actuator dynamics 
to the structural dynamic models. This will allow the mass and stiffness as 
well as the damping rate of the damper to be design variables. Thus this will 
be the actuator design phase. 

Next, a 2-dimensional structural model which has a higher modal density 
will be developed. 



NONCOLOCATED SENSORS AND ACTUATORS 
ADDITION OF ACTUATOR DYNAMICS 
2-DIMENSIONAL STRUCTURAL MODEL 

Figure 7 



12 



REDUCED ORDER ADAPTIVE REGULATION 
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REASONS FOR ADAPTIVE REGULATION 

Uncertainty in modal frequency (and mode shape) 

prespecification 

Inability to perform identification prior to 

deployment 

Time-variations in system dynamics during operation 

(especially construction) 

NASA BEAM MODE FREQUENCIES 
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REASONS FOR MODAL DESCRIPTION 

Decoupled dynamics convert multi-input, 

multi-output problem to several 

single-input, single-output problems 
Parallel computation for real-time control implementation 
Common form of models of large flexible structures 
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MODAL DESCRIPTION OF FREE-FREE BEAM 
y (x,t) = ^ y t (x,t) 

i l 

yj(x,t) = 4>j (x) vj (t) 

"i (t) + w iS (t) = Fi (t) 
Fj (t) = zV(xj) f.( X ,; t) 



j 



1 W Y V -J 



»l (x) * [cosh WjX + cos.-Wi-x 

(cosh Wjl - cos Wjl) 



, ., ,,, . lx (sinh w<x + sin w^)] 
(sinh W| 1 - sin Wjl) V l 



Note: +j dependent on Wj 



y : vertical deflection 

x : horizontal beam point location 

t : time 

yj : ith mode deflection 

*i :' ith mode shape (characteristic function) 

*! : ith mode amplitude 

W| : ith mode natural frequency 

F| : ith mode forcing function 

fj(Xj, ■): Jth point actuator force applied at x. 

1 : beam length 
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QUESTIONS 

• How do currently available single - input, single-output 
adaptive control schemes behave when applied as modal 
controllers of flexible structures? 

• Can adaptive regulation (stabilization) be achieved with 
reasonable control levels despite inexact mode frequency 
(and shape) prespecification? 

• How close must initial estimates be? 

• Is adaptive regulation "better" (i.e. faster, closer to 
some desired behavior, ...) than fixed or gain scheduled 
control of comparable complexity designed for some bounded 
uncertainty in mode frequencies(and shapes)? 

• What (if any) augmentations to currently available 
adaptive controllers improve their performance with 
application to flexible structures? 
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ADAPTIVE MODAL REGULATION CANDIDATE 

Given: expansion basis Uj}, finite expansion limit N, sensor 
and actuator locations, and modal control objectives 

(1) Apply previously calculated actuation forces f,(x, , t) 
and sense y(x, t) 

(2) Process sensor data to estimate M .modal amplitudes 
yj(x ,t) 

(3) Select C. modes requiring control 

(M) Process applied forces fj(x.,t) to determine achieved 
modal forces Fj(t) 

(5) Improve the identification of the discretization of the 
^ dynamics (i.e., essentially the modal frequencies w t ) 

(6) Parameterize modal controllers using; current plant 
parameter estimates, to meet modal performance objectives 
and calculate F^t + 1) 

(7) Convert desired modal forces Fj(t + 1) to actuator 
commands f j (Xj ^ t + 1) 

(8) Repeat (1) - (7) 

Source: Johnson, "Adaptive modal control of large flexible 

spacecraft, " J . of Guid. a nd. Control, July-August 1980, 
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PROBLEMS 

Finite number of modes accurately describing system N > 
number of modes identifiable M due to number of point 
sensors > number of modes controllable C, i.e., the reduced- 
order adaptive identification/control problem. 

Disregard of coupling due to inexact modal shape +j 
specification. 

Unbounded control effort request due to momentary uncontrol- 
lability of estimated plant parameterization. 

Lack of, guaranteed .parameter 'ideritif lability without 
sufficient excitation. 

Misinterpretation of deflection due to short-term distur- 
bances as caused by parameter estimate inaccuracy. 

Meaningful prespecification of modal control objectives 
(e.g., pole placement for Quadratic cost function 
minimization) given modal frequency and shape uncertainty. 
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ALTERNATE ADAPTIVE STRATEGIES 

Asymptotic feedback matrix synthesis can bypass momentary 
estimated plant uncontrol lability. 
Source: Kreisselmeier, "Adaptive control via adaptive 
observation and asymptotic feedback mat irx synthesis/' 
IEEE Trans, on Auto. Control , August 1980. 

Avoid modal decomposition and treat as multi-input 

(C actuators), multi-output (M sensors) problem to 

accomodate "modal" coupling with loss of computational 

advantages. 

Source: Goodwin, Ramadge, and Caines, "Discrete-time 

multivariate adaptive control, " IEEE Trans, on Auto . " 

Control , June 1980. 

Use direct adaptive controllers to bypass perturbation 
necessity for identif lability. 

Source: Johnson and Tse, "Adaptive implementation of one- 
step-ahead optimal control via input matching, " IEEE Trans . 
on Auto. Control , October 1978. 

Use off-line parameter estimation with verification prior to 
alteration of gain schedule. to avoid reaction to non- 
parameter-estimate errors such as deterministic disturbances. 
Source: Hall, "A learning control system extension to 
the modal control of large flexible rotating spacecraft," 
Proc. 1979 AIAA Guid, and Control Conf. , August 1979. 

Limit controller parameterizations to positive real operators 
with subsequent loss in objective range in order to assure 
stabilization despite reduced-order usage. 
Source: Benhabib, Iwens, and Jackson, "Active vibration 
control of a flat plate using model, reference adaptive 
techn i ques , " Proc. 2nd VPI&SU/ AIAA Symp. on Dyn, and Control 
of Large F lexible Spacecraft, June 1979, 
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SPECIFIC QUESTIONS TQ BE ADDRESSED 

How many modes (with N ■'= M = C) are required in a modal, 
simultaneous identifier-controller, initialized with 
the ideal description of a free-free beam, to satisfactorily 
regulate Langley's flexible beam? 

Using only displacement measurements what order actuator- 
input, sensor-output, mat rix-ARMA model is satisfactory 
for nori-modah -multi variable adaptive regulation of 
Langley's flexible beam? 

By describing the effects of reduced-order modeling as 
a perturbation to" the time-varying, nonlinear, parameter- 
estimate-error system, can the theoretical limitations of 
reduced-order adaptive control, especially in terms of 
initial estimate accuracy, be interpreted? 



Using a singular perturbation separation of the modes of 
a flexible structure, can the effects of reduced-order 
adaptive control be quantified? 
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MILESTONES 

1. Select and prepare a mission model 

2. Write a set of equations of motion 

3. Determine elastic constants of model components 

4. Write a computer program for open loop input 

5. Set up control block diagram 

6. Simulate the dynamic system on the computer 



MISSION MODEL 

The Orbiter and the SEP solar array are represented by the model shown in figure 
1. The coordinate system is shown on the diagram. 



~2L 



zz—1 



— =R7^ 




Astromast 



Solar panels 



-&r-7 




Shuttle 




(a) Real system. 



(b) Ideal system. 



Figure 1 
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KINEMATICS OF RIGID BODIES 



V„ = V + V , 
B A B/R 

V B/R = w x ^ 



W = rotation vector 



r = position vector 
*B = ft ( V 



SEP SOLAR. ARRAY MODEL 
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MAST MODEL 

Assumptions: 

1. The total mass of the mast is concentrated at each level. 

2. Lateral members are relatively much stiff er than longitudinal ones. 

3. The deformations in the mast do not depend on axial forces present in 
the londitudinal members. 
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m, • 1 



J"i 



Asiromciit 



77)7777 

Ho6c\ 

:°+ 

Mast 



id ealned 
Mast 



hi} = tJilidtt + d * uJ 

2BI HI I** 

neXIBILlTY MATRIX » HI 
STIFFNES5 MATRIX « CK3 = CO 
/MASS M/ITKiX CD/AGONAL), f/vj] 
fOKCH MATRIX «. CFJ 

[MXK-) + [K]£>0 « o 

Figure 3 
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MODAL ANALYSIS OF ASTROMAST 
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Figure 4 
MODAL FREQUENCIES 
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INTRODUCTION 

Complete decoupled control requires the number of control actuators to equal the 
number of modes in the model, which is a basic limitation in applying decoupling 
theory to the control of large space structures. Complete decoupled control is 
usually not achievable in practical application because a large space structure may 
have an infinite number of flexible modes; hence, procedures must be developed which 
maintain control of the structure with a small number of control actuators. The 
present analysis presents techniques which use decoupling theory and state-variable 
feedback to control the pitch attitude and the flexible-mode amplitudes of the beam. 
Approximations are incorporated into the decoupling procedure to permit control with 
a small number of actuators. The approximations involve adjustments in the control- 
influence coefficients and in the feedback gains which produce simplified procedures 
for achieving overall control of the system. 
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EQUATIONS OP MOTION 



Figure 1 shows the linearized equations of motion used for a decoupled-control 
analysis of a 100-m long thin, flexible beam in low Earth orbit. The equations are 
in modal form where n represents the number of flexible modes included in the model 
Although not required, the damping term 2C0)A n is included in the model. The top 
equation represents the rigid-body (pitch) mode and includes the gravity-gradient 
effect, where co c is the orbital frequency. 
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dt 
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Figure 1 
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DECOUPLED CONTROL EQUATIONS 

Figure 2 shows the equations in state-vector form where the states x are the 
modal amplitudes and rates. The output equation y = Cx represents the states to be 
decoupled. In the decoupling control law v is the input command vector and F 
and G are the feedback and feedforward gain matrices, respectively. The output is 
related to the input through the transfer function H(s). The decoupling procedure 
determines the F and G matrices in a manner. such that the transfer function is 
diagonal and nonsingular, thus providing independent control for each of the 
decoupled (output) variables. 

SYSTEM XCAjB.C) 

x = Ax + Bf 9 g s:- Cx 

DECOUPLING CONTROL LAW 
J ss Fx 4- Gv 

LAPLACE TRANSFORM SOLUTION 

^(s) ~ Hcs)vC-s) 

Where W(s) ~ C(sl-A - BF)~ ! BG 
is the wxm transfer -Function, 

For decoupling , H(s) is diagonal and 
nonsingular. 

Figure 2 
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EXAMPLE EQUATIONS 

Figure 3 presents example equations for a four-mode model (pitch plus three 
flexible modes) with four control actuators. The numbers in the A matrix are the 
frequencies squared of the various modes. The damping terms have been omitted 
because they are not required in the solution for decoupled control- The control- 
influence matrix B depends on the location of the control actuators and the modal 
shape functions . 
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DERIVATION OF DECOUPLING CONTROL LAW 

Figure 4 is an example of a derivation of the decoupling control law for the 
model of the previous figure. Inasmuch as the number of control actuators equals 
the number of modes, complete decoupling is obtained. That is, each of the modes can 
be independently controlled. Obviously, if the number of actuators is reduced, some 
of the modes will remain uncontrolled. These uncontrolled modes will be influenced 
by effects because of control of the decoupled modes. Subsequent results will 
illustrate several methods for reducing or eliminating these undesired effects. 
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EQUATIONS IN MATRIX FORM 

Figure 5 shows the equations for the control law in matrix form. Note that the 
feedforward matrix G is merely the inverse of the control-influence matrix. Each 
column of G can be multiplied by a constant to change the output sensitivity. Note 
that the matrix M is composed of the modal frequencies to and the damping 
ratios £. These are closed-loop quantities and their values can be selected to 
provide feedback gains for desired dynamic-response characteristics. Note that the 
feedback gain matrix can be separated into two parts, that part F" which deletes 
the unaugmented dynamics and F" which incorporates the dynamics selected for the 
closed- loop system. 

v - Mx + Mf 

or ? = -NT'Mx + fSf'v 
.Providing N is not singular^ then 

where 

F= -N~ l M = F'+ F" 
G - N"' 

are the required feedback and feed-forward 
matrices for the decoupling Control law. 

Figure 5 
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PITCH COMMAND FOR COMPLETELY DECOUPLED MODEL 

Figure 6 is an example of an instantaneous pitch command to change the pitch 
attitude by 0.01 rad. The quantities A,, Aj, and A^ are the modal amplitudes; 
fj_ and f2 are the actuator forces in newtons . As shown, the system is completely 
controlled in that the flexible modes are not affected. The results apply to the 
condition that complete feedback of all modal amplitude magnitudes and rates is 
required. Also, the effect of unmodeled modes, not included in the computer simula- 
tion, is not considered. 
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Figure 6 
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ZERO COMMAND FOR COMPLETELY DECOUPLED MODEL 

Figure 7 is an example of an instantaneous zero command to null initial dis- 
turbances of 0.01 rad in pitch and 0.01 m in the modal amplitudes. The values 
selected for the closed-loop dynamics of each mode were 0) = 0.1 rad/sec and 
Z, = 0.5. This gives similar responses in all modes leading to equilibrium in 
about 60 sec. 
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ONE CONTROL ELIMINATED 

Figure 8(a) is an example of a pitch command where the f 3 control (actuator) 
was eliminated from the completely decoupled system. The same F gains derived for 
the completely-decoupled model were used, except that the gains for the f 3 control 
were deleted. As shown, this condition leads to uncontrolled first and third 
flexible modes which are caused by oscillations in the f^ and f 4 actuators. 



Figure 8(b) shows how overall control of the previous case can be obtained 'by a 
simple gain adjustment. The example is for a zero command and shows that the initial 
disturbances are controlled (nulled) after about 100 sec. The gain adjustment 



required in this case was a change in F 
the feedback matrix) by a factor of ten. 



1,8 



(element in first row, eighth column of 



More than one control can be eliminated from the completely decoupled system; 
however, as more controls are eliminated, a larger number of adjustments in the feed- 
back gain matrix are required (involving many trial and error attempts) . 
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ZERO COMMAND WITH TWO ACTUATORS 



fi -^ 9Ure / presents an example of a zero command where the pitch mode and first 
flexible mode are decoupled with two control actuators. This leads to incomplete 
decoupling due to the reduced numbers of controls. However, as shown, except for 
some initial effects m the modal responses, overall control of the system is main- 
T^u \ 1S a Special case inasmuch as the actuators are exactly at each end 

of the beam. Gain information derived from this special case can be useful for 
adjusting gains for other control arrangements which lead to uncontrolled modes 
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CONTROL-FORCE REQUIREMENTS 

Figure 10 illustrates the control-force requirements for the two-control model. 
The values shown apply to a pitch command of 0.10 rad. It is important to note the 
high dependence of control force on the value selected for the closed- loop pitch 
frequency. At the larger values of OJg the forces become prohibitive. The forces 
for the higher-mode models are somewhat lower because the F matrices require 
higher-order dynamics to calculate the gain values for these models. 
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EFFECT OF PARAMETER UNCERTAINTY 

Figure 11 shows the negligible effect of parameter uncertainty on a pitch com- 
mand; the solid curves represent the no-error case. The dashed curves represent a 
typical result from a number of cases where random errors of ±10 and ±20 percent 
were incorporated into the control-influence matrix B. For these cases, the feed- 
back gains calculated for the original B matrix were used. The uncertainty in 
beam model parameters was also found to have no appreciable effect on the decoupling 
process. These uncertainties included errors of 20 percent in the modal frequencies 
and large changes in the modal damping factors (up to 500 percent of the nominal 
values) . 
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PITCH COMMAND WITH TWO ACTUATORS OFFSET 

Figure 12 shows an example of a pitch command where the pitch mode and first 
flexible mode are decoupled with two control actuators. This incomplete-decoupling 
case differs from the previous case in that one of the actuators is not exactly at 
the end of the beam. As shown for this offset case (which is the general case) , the 
required decoupling is achieved for G and A]_; however, the other two flexible 
modes remain uncontrolled. This result is described in the next figure. 
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FEEDBACK GAINS REQUIRED TO CHANGE DYNAMICS 

Figure 13 shows the feedback gain matrices derived for the two cases, one for 
offset actuators and one for actuators exactly at each end of the beam. Only the 
first four columns are shown for F" ,, inasmuch as the last four repeat by 
replacing S with R and Q with P. As shown for this matrix, fourth-order 
dynamics are required to calculate the gains. The constants P, Q, R, and S 
depend on the closed-loop values selected for oj and £,. The zero columns in 
F offset' which provide no control for the second and third flexible modes, result 

from B-matrix columns of unequal magnitudes. The f" matrix is full order in 

exact 
that no zero columns exist. This condition occurs only when the values, in the col- 
umns of the B-matrix have exactly the same magnitude. The ratios between the 
columns in F gxact can be used to adjust the feedback gains for the zero columns 

of the offset case. In applying the ratio method, the third column of f" is 

offset 
adjusted so that its ratio with the first column in F ff is the same as the 

respective columns in F g Xact - The same procedure is used with respect to the 

fourth and second columns, the seventh and fifth columns, and the eighth and sixth 

columns. A result of this gain-adjustment procedure is shown in the next figure. 



'exact ~ 



84oo ,250 
•8*oo 25o 

250 ZSO 
-25o Z5o 



S a - Q e <«W 



5, -Q,fa> lu , 





o 



o 
o 



o 
o 



o 
o 



o 



W'a — tO. 



= s 4 + p $ 3 + QS Z + R S + S 



F offset = 



84oo 25Q 
■84oo ZSo 



E©o oo D©ooo 
O E, o o o D, O o 



% ^ cZ. 



sHi^wstw = s*+D*+-E 



Figure 13 



41 



ZERO COMMAND WITH TWO ACTUATORS OFFSET 

(FEEDBACK ADJUSTED BY RATIO METHOD) 

Figure 14 shows an example of a zero command where the feedback gains were 
adjusted by the ratio method previously described. As shown, the adjusted gains 
provide sufficient control for all modes. There are small initial effects, but they 
damp out after about 200 sec. 
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UNCONTROLLED MOTIONS FOR FOUR-MODE MODEL 

Figures 15(a) and 15 (b) illustrate the: use of Kalman regulator theory in esti- 
mating the modal amplitudes. Similar results are obtained for the amplitude rates. 
These estimates are required for the feedback control inasmuch as the actual modal 
values cannot be measured directly. The first plot shows the uncontrolled motions 
resulting from, small initial distrubances . The second. plot shows how well these 
motions can be estimated with measurements from one pitch-attitude (star tracker) 
sensor located at a point one quarter, the distance, from the end of the beam. This 
plot shows that perfect estimates are established at a time of about 60 sec, where: 
the magnitudes of the oscillations are equal to those of the first plot. (Oscilla- 
tions are actually present in the estimates of A lf but are hot discernible due to 
the large scale,) In. this one-sensor case, if the control actuators are turned on 
before good ..estimates are obtained, excessive control forces would be required. 
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ZERO COMMAND WITH THREE SENSORS 

'Figures 16(a) to 16(c) show zero-command results for a three-sensor case. The 
sensors were located at the middle and quarter points of the beam and the control 
actuators were at each end of the beam. In this case, the initial estimates were 
adequate so that the actuators could be turned on at t = without requiring 
excessive forces, as seen in the last plot. These forces were about six times those 
required for the case with perfect initial estimates. With only two sensors 
employed, the forces increased substantially (to about 65 N) , which would require 
turning off the actuators until good estimates are produced. 



These Kalman- regulator results apply to the case where unmodeled modes are not 
considered in the computer simulations. In any practical application, unmodeled 
modes cause observation and control spillover effects which can lead to instability. 
Additional analysis is required to include these effects in determining the overall 
stability and control of the system. 
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1 . Introduction 

The feasibility of constructing large space structures such as the solar power 
system, the SEP (solar electric power) array, or numerous other space structures 
currently under study will undoubtedly be realized in the future. The existence of 
such structures will depend to a great extent on the ability to analyze and control 
large scale structures. The dynamic behavior will require algorithms that are effi- 
cient for handling thousands of differential equations. These types of problems ■ have 
been included in disucssions on weather prediction by Bellman where tens and hundreds 
of thousand partial differential equations are required for local weather phenomena. 

Currently available algorithms will most likely be unable to handle large scale 
systems because of storage requirements, execution time, or numerical accuracy. The 
investigation of algorithms for large scale problems is part of the research being 
carried out by this investigation. The primary objective of the research is to 
obtain an efficient algorithm for decoupling a large set of differential equations 
into subsets of differential equations that can be numerically integrated. The 
solution of the large set is then constructed from the subset solutions. 

Problems Under Consideration 



1. Spectral factorization 

2. Decoupling of a large number of differential equations 

3. Eigenvalue-eigenvector subroutines 

4. Matrix polynomials 

5. System identification 

2. Sign Algorithm 

The sign algorithm was first introduced by Roberts as an algorithm to compute 
the solution to the algebraic matrix Riccati equation given in equation (2.1) . It 
is not difficult to show that the matrix given in equation (2.2) can always be 
associated with the matrix Riccati equation. The two solutions to equation (2.1) 
are given in equation (2.3) where $. . are partitioned blocks of the eigenvector 
matrix of A as given in equation (2^4) . The sign algorithm given by Roberts is a 
Newton procedure where S. in equation (2.5) will converge to the sign of A. 

Let A be 2nx2n with n eigenvalues having Re (X. ) > and n with 
Re(X i ) < 0. If J is as defined in equation (2.6) then sign (A) is as given in 
equation (2.7). It can be shown that the inverse of S + J contains the two solu- 
tions to the algebraic matrix Riccati equation (see eq. (2.8)) and is related to the 
eigenvector matrix of A as given in equation (2.9) . 

The T matrix will block diagonalize A when applied to A as a similarity 
transformation. The resultant of the operation is shown in equation (2.10) . The 
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spectrum of 



will contain only eigenvalues with 



_, -, Re (A.) > 0, whereas that of 

will have Re(A-) < 0. The transformation on A is a spectral decomposition 



,2 
and is useful in solving a system of differential equations. Let x(t) be a 2nxl 

vector satisfying the state equation (2.11). If z(t) is as defined in equation 

(2.12), then equations (2.13) and (2.14) hold where z,(t) and z 2 (t) are nxl 

vectors which are partitions of the 2nxl z(t) vector. The system response vector 

x(t) can be recovered by using equation (2.15) or (2.16). The procedure can 

obviously be extended to more than two diagonal blocks. 



Knowledge of the sign matrix of A is sufficient to define a set of eigen- 
projectors P^ and P 2 for a simple spectral decomposition. If P^ and P 2 are 
as defined in equation (2.17), then A-^ and A«» as given in equation (2.18), will 
have spectrums p(A 1 ) >^ and p(A„) 
also decompose x(t) into two 2nxl 



<_ 0. The eigenprojectors 
vectors x + (t) and x~(t). 



P, and P 2 will 
The vector x + (t) 



will be a function of exp(A.t) 
modes with Re(A^) <_ 0. 



with A. having 



Re(A ± ) > 



and x~(t) will contain 



The eigenprojectors introduced in the previous paragraph were computed from the 
sign function computed on the basis of splitting the spectrum along the jco axis. 
The bilinear transformation of equation (2.20) will permit a spectral decomposition 
with respect to circles of radius p chosen as to split the spectrum by eigenvalue 
magnitudes. The previous analyses still hold and, as before, several transformations 
may be used to decompose the spectrum into more than two domains. 

The above decoupling procedure has several disadvantages. The first of these 
is the requirement for inverting 2nx2n matrices which is costly and inaccurate when 
n is large. In addition, several large matrices must be stored. An alternate 
approach is to use a standard eigenvalue-eigenvector routine to compute S. There is 
very little to gain in this approach, since the time domain solution can then be 
determined without the use of the decoupling. 
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(2.5) 



.w^^v*^ 



S = A. 



(2.6) 



J = diag[l - i] 



I-..is.,nxn .identify 
matrix 



(2.7) 



S = $J$ 



-1 



S is 2hx2n 



(2.8) 



T = (S + J) 



-1 



I - . -E. 



T is 2nx2n 
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T = 
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x(t) = Ax(t) + Bu(t) 
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(2.14) 



z 1 (t) = Ag^^t)' + B^it) 



z 2 (t) " A B2 Z 2 (t) + B 2 u(t) 



(2.15) 



x(t) = T z(t) 



(2.16) x x (t) = ^E^t) - R z 2 (t)] 



x 2 (t) = |[R 2 1 z ;L (t) - z 2 (t)] 



(2.17) 



P x = | (I + S) 



P„ = - (I - S) 
2 2 



P. + P 2 = I 



(2.18) 



APj^ = A r AP 2 -= A 2 -.-.A ■'.+ A 2 =• A 



(2.19) 



P 1 x(t) = x + (t) P ? x(t) =x"(t) ; X + (t) + x~(t) = x(t) 



(2.20) 



A„ = (A - pi) (A + pi)" 1 .. ,p < I A, 
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3. Matrix Polynomials 

The difficulties encountered in the previous section are less severe if the A 
matrix is in or can be transformed into the block companion form and certain algorithms 
can be developed. As an example, assume that A has the form given in equation (2.2) 
which can be transformed into the second order companion form of equation (3.1) by 
a Krylov transformation given in equation (3.2). A typical Krylov matrix is given 
in equation (3.3) where it is assumed that A l2 is invertible. The matrix polynomial 
of equation (3.4) will have 2n latent roots A. and 2n latent vectors y^ where 
the latent roots A. are the eigenvalues of A and A . The latent vectors y. 
are the 2n vectors of ^ and $_ where $ is the eigenvector matrix of A . 
Since A(A) is of the form of the differential equation for a finite element model, 
the second order matrix polynomial is usually of interest in structural problems. 
The coefficients of the matrix polynomial, the latent roots, and the latent vectors 
of A(A) require less storage than that required for storing A c , the eigenvalues 
and the eigenvectors. This approach is therefore of interest because of storage 
requirements. The disadvantages of formulating the problem in terms of matrix poly- 
nomials is that efficient algorithms for computing the latent roots and the latent 
vectors of A(A) do not exist. This is an area that is currently under investiga- 
tion . 

A substantial decrease in storage requirements can be obtained by transforming 
A (A) into a higher order polynomial with matrix coefficients of lower order. Assume 
that A(A) is. nxn and that A is 2nx2n. If 2n is divisible by m in an integer 
sense, then an mth order polynomial can be found when A(X) has distinct latent 
roots. The matrix coefficients of A (A) , as given in equation (3.6) , will now be 
(2n/m) x (2n/m) matrices. The latent roots of A(A) and the latent vectors are 
sufficient to define the eigenvalues and eigenvectors of A where A is 2nx2n as 
given in equation (3.7). 

The computation of the eigenvalues and eigenvectors of A or the latent roots 
and latent vectors of A (A) is not necessary if the time domain solution to the 
system differential equation is sought. 

The Krylov type transformation required to construct the higher order poly- 
nomial is also used to transform x(t) to a new vector z(t) , as given in equation 
(3.10), where z(t) will be partitioned as in equation (3.11). It can be shown 
that z 1 (A) will be defined by equations (3.12) and (3.13) with zj(t) given by 
equation (3.14). x(t) can then be constructed from the Zi(t) with x(t) = K z(t) . 

A second method of analysis is to use block matrices throughout. Let <j). be 
solutions to the mth order polynomial, then A(A) can be factored as in equation 
(3.15) . Thus, z.. (t) can be found from equations (3.16) and (3.17) . The z. (t) 
required for constructing x(t) are found in a similar manner with the proper 
modifications in the matrix coefficients of equation (3.17) . 
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The procedure given has several associated problems which must be solved prior 
to implementing the above algorithms. The first of these is that there does not 
exist an efficient algorithm for the Krylov transformation to construct the mth order 
polynomial from a lower order polynomial. Secondly, there does not exist at this 
time an efficient algorithm to compute the latent roots of A (A) or the latent 
vectors. The factorization procedure expressed in equation (3.15) is also an area 
for development of an efficient algorithm. The development of these algorithm is 
presently under study. 
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(3. 9} [x(t)3 T = [x T (t) x T (t)3 

(3.10) z(t) =' K x(t) '"'■' : '..'.-' V 

(3.11) z T (t) ■■= [z x (t) -z 2 (t) ■■■: . . z m (t7| z i (t) is 2n/mxi 

(3.12) A(X) z (X) = B(X) U(X) MX) is 2n/mX2n/m 



2n P . n 



(3.13) z (X) = [a(X)]~" l B(X) u(X) = S. (A _ X ) B( V U(X) 

i=l i 



2n ^ t 

(3.14) z (t) = S P B(X ± ) [ exp [X (t - T)] u(T) dx 
i=l JO 



(3.15) A(X) = (IX - Q 1 ) (IX - Q 2 ) . . . (IX - Q m ) 



m B((j).) m 
(3.16) [a(X)]"b(X) = S — - = S P iQ 

ill <*± - v 



m 
(3.17) z (t) = S P iQ exp iQ^t - T)]u(T)dT 
i=l 



4. System Identification 

Numerous algorithms have been written for the identification of a system and for 
parameter estimation. The method of least squares solution of a set of overdetermined 
equations appears to be the most widely used algorithm, although there are other 
methods that are acceptable. The least square algorithm basicly takes the discrete 
time data from the system response with a known input to define an overdetermined set 
of equations. As an example, free response, is given in equation (4.1) for a time 
invariant system having a typical set of equations for a system with no input. The., 
matrix is the state transition matrix which will be invariant for all t. Equa- 
tions (4.2 to (4.5) can then be established. 

The unknown parameters in are then computed by a least squares algorithm 
such as given in equations (4.6) and (4.7). This approach is riot recommended but is 
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given only to illustrate the simplest procedure. The singular value decomposition 

would be preferred for determining the vectors 0. . 

1 

The complete set of row vectors 0. would be identified and would then 
be known. Since is discrete, the slate transition matrix will have the form 
given in equation (4.8). The eigenvalues and eigenvectors of could then be 
identified. It has been assumed that all computations were made with no error 
although this is an ideal situation. It must also be assumed that the sampling 
frequency of the system response has been sufficiently high to satisfy the sampling 
theorem. 

The numerical problems encountered in the decoupling algorithm will be present 
in the identification procedure. The storage ; requirements of the algorithm dictate 
that a large machine be available, and the algorithm will probably require multiple 
precision. The execution time for a system with 1000 eigenvalues on a typical 
machine may be an hour. 

The research has concentrated on the decoupling algorithm, eigenvalue and eigen- 
vector algorithms, matrix polynomials, and identification procedures. The decoupling 
algorithm as presently used requires inversion of large matrices which is inefficient. 
Experience with applying the algorithm to a free-free beam indicates that numerical 
errors are severe when the beam has several hundred nodes due to the condition number 
of the matrix or the range of the elements in the matrix. 

Numerical errors have also been apparent in the use of the eigenvalue-eigenvector 
subroutines from EISPACK on relatively small matrices, 30 x 30. The free-free beam 
problem with twenty-one nodes has been decoupled with the algorithm discussed in the 
next section, and it has been necessary to use an eigenvalue-eigenvector subroutine 
with all computations in double precision on a 30-bit machine . 

Research on matrix polynomials has been underway, and some of the results were 
presented in Section 3 of this report. This mathematical development of algorithms 
for analyzing matrix polynomials has not received the attention that seems to be 
warranted. Some aspects of the use by matrix polynomials will be covered. 

System identification of large scale systems has not received much attention in 
the past, since most practical systems are of order 100 or less. Identification of a 
system with 200 or more parameters is not a commonly occurring computational task. 

The use of a matrix polynomial in identification has received very little, if 
any, attention in the past. The identification algorithm by quadrature methods uses 
scalar functions for determination of the parameters. The question arises as to 
whether matrix functions can be utilized rather than the scalar functions. This 
question and many other remain unanswered. 

(4.1) x(k + 1) = 0x(k) T = sampling time: k = t/T 
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(4.2) 



Xl (i) =e llXl (o) + e i2 x 2 ( 0) + ... +9 ln K n (o) 

x l( 2) = 11 x 1 (l) + 12 x 2 (D + .-• + ln x n (1) 



x l (m) = ^lV™" 15 + 12 X 2 (m_1) + '" + ln X n (m_1) 
(4.3) x T (k) ■= [x 1 (k) x 2 (k) ... x n ( k )] T 



(4.4) 



x 1 (l) 
x 1 (2) 



x (m) 



x 1 (0) 
x (1) 



x 2 (0) 
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x (m-1) x (m-1) 



x (0) 
n 

x (1) 
n 



x (m-1) 



9, 
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In 



(4.5) e 
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n 



(4.6) Y = r X 



(4.7) 0t = (R T R) _1 R T Y 



(4.8) = $e 



At*-1 



At r 1 2 n -i 
(4.9) e = diag Le e . . . e J 
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INTRODUCTION 

Controlling large, flexible objects in space is currently a "hot" topic in the 
research community. Many "solutions" have been proposed, most of which have 
either not worked in practice or have flaws which come to light when one considers 
implementation. At Honeywell, we have done extensive research on implementable 
control designs, which has served to guide the research to be performed. One of 
the critical problems we see is in obtaining and then properly using accurate 
information about the structure to be controlled. 

By way of motivation, I will first give a brief review of our view of LSS control, 
describing the problem and what's new. Following this, the results of applying a 
new identification measure to a sample LSS will be presented. These results indicate 
the promise/problems of on-orbit identification. I will conclude by summarizing the 
line of research we are pursuing. 

8 Motivation 

e Identification 

Our Program 

e Expected Results 

Figure 1 
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THE CONVENTIONAL PROBLEM 

In the past, spacecraft flexibility has been handled by requiring the flexible modes 
of the spacecraft to occur well outside the bandwidth of the control system. Where 
this did not occur naturally, system modifications were made to ensure that it did 
happen. By then "rolling-off " the control system, interactions with the flexible 
dynamics were avoided. 



TRANSFER FUNCTION 
MAGNITUDE 



STRUCTURAL MODES 
(FLEXIBILITY EFFECTS, ETC.) 



CONTROLLER BANDWIDTH 
(POINTING, SHAPE REQUIREMENTS, ETC.) 




• • • 



»-CJ 



SOLUTION: "ROLL OFF" CONTROLLER TO AVOID EXCITING HIGH FREQUENCY MODES. 

Figure 2 
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THE PROJECTED PROBLEM 

As spacecraft grow and become more flexible and as control bandwidth increases due 
to increasingly stringent control requirements, the solution of forced separation 
is becoming intolerable. Talk is heard of hundreds of flexible modes within the 
bandwidth. This overlap characterizes what we call "the LSS control problem." 

The solution is then clear: one must "actively control" the modes in the pass band 
of the controller and roll-off (while maintaining stability) the higher frequency 
modes. 



STRUCTURAL MOOES 




*-cj 



CONTROLLER BANDWIDTH 



SOLUTION: 1) ACTIVE STRUCTURAL CONTROL OF MODES 
WITHIN CONTROLLER BANDWIDTH. 

2) "ROLL OFF" CONTROLLERS TO AVOID EXCITING 

STRUCTURAL MODES BEYOND CONTROLLER BANDWIDTH. 

Figure 3 
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INFORMATION TRADE-OFF 

There are two basic approaches to achieving these. The first is a robust controller 
based on colocated actuators and sensors using rate feedback. This approach is 
extremely insensitive to model information requiring only crude mode shape and 
frequency information and only a limit to damping. However, the amount of modal 
damping achievable this way is limited and relatively high bandwidth actuators, 
sensors, etc., are required. 

In contrast, many methods have been proposed which achieve better damping and require 
less bandwidth from components. However, invariably, these require accurate model 
data. 

• Robust Control 

Insensitive to Model Information 

- Achieved Damping Limited 

- Requires High Bandwidth Intelligence 



High Performance Control 

- Achieves High Damping 
Requires Less Bandwidth 

- Very Sensitive to Model Data 

Figure 4 
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WHY IDENTIFY 

In general, pre-*launch data is not accurate enough to permit design of high 
performance controllers. Thus, if the benefits of these are to be realized, system 
identification must be performed. 

• Want ienefits of High Performance Control 

• Pre-launch Data Poor 
Only Two Options 

Adaptive 
- Identify and Reconfigure 
§ Either requires identification 



Figure 5 

IDENTIFICATION 

Two basic methods of system identification are currently popular in the aerospace 
community: those based on Maximum Likelihood Estimation (MLE) and those based on 
the frequency domain (FFT) . Some recent advances in the MLE analysis are very useful 
for our type of problem. We have applied these to a sample problem, the results of 
which I would like to present. 

• Maximum Likelihood Estimation (MLE) Based 



• Frequency Domain Based 

Figure 6 
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MLE 



The MLE method is based on finding which model from a (possibly infinite) set of 
models comes closest to the . true system. Closest, for the MLE method, is defined as 
minimizing the negative log likelihood function, which can be computed from the 
Kalman filter based on a model M . 

a 

Of course, LSS are infinite dimensional and therefore no finite dimensional model can 
ever perfectly match the true system. This mismatch has caused problems in the past 
and until recently no method of analysis was available. 

x(k+1) - AWx(t) + B.u(t) + L.£(t) 
Z(t) - C»x(t) ♦e(T) 

Truth: M. - ( A., B., C, U H., 9.) 

Model: « a - JA a , B a ; C a , L a , Ea, 9 a | 

6oal: Find model M. "closest* to FW. 

Closest - Minimize the negative log likelihood function 
(computed from KBF for M a ) 

Figure 7 
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DISTANCE MEASURE 

Based on recent work by Baram, extended by Yared (Ph.D. thesis, MIT, 1979), a 
distance measure relating the information favoring one model over another can now be 
computed. Using this, the distance from a model to the true system can be computed. 

Note that this assumes the true system to be known. It is, however, useful for 
analysis purposes for predicting the behavior of the MLE method under various model 
order, initial condition, and noise conditions. 

Io ■ Expected value of the conditional likelihood function 

WITH RESPECT TO THE TRUE SYSTEM. 



K°l,«2) * la 2 " 1«1 " DISTANCE MEASURE 

Information favoring model a x over a 2 . 
KV) ■ Distance from true system to model «. 

Figure 8 



64 



SAMPLE TETRAHEDRAL TRUSS STRUCTURE 

We have applied this technique to a sample structure used in the ACOSS program known 
as the tetrahedral truss. Test signals (forces) and displacement measurements were 
assumed in one leg of the truss. Noise sources were included as shown. 

Tests were then conducted to determine the distance between the assumed true system 
and a model. For all tests, the model used for identification contained a single 
mode while the number of modes used to represent the true system was either one or 
nine. In all tests, only one parameter was assumed unknown. 



NODE NUMBER 




MEMBER 
NUMBER 



Forct Actuator 
Displacement Stnsor 



Notes: 



Plant Disturbance 
Sensor Noise 
Test Input 
Damping Ratio 



.005 rms (produces .001 at sensor) 

.001 rms 

.05 rms (produces .01 at sensor) 

,5X 



Figure 9 
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DISTANCE FROM, ONE MODE TRUTH FOR 0) ID 

In the first experiments, the objective was to identify the frequency of a model 
assuming, the damping to be known. • If . the true system and identification model match 
in order, the distance is as below. 

Note that the distance is zero at the : true -value of the parameter. Also, the valley 
is very steep, indicating good local convergence but requiring good initial guesses. 



e.2G 



8,15 



0.10 



0.05- 



8.00 




1.5 2.0 2.5 3.0 

Estimtt of w. 

True Value 

Figure 10 



4.0 
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DISTANCE FROM 9 MODE TRUTH FOR CO ID 

If the true system is larger, the distance becomes large and no value of co reduces 
the distance to zero. In addition, multiple minima occur. 

48.5 



49. e 



39.5 



39.8 



38.5 




1.8 



- 



True Value 



2.0 2.5 

EstlMte of Wj 

Figure 11 



3.0 



3.5 



4.8 
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DISTANCE FROM ONE MODE TRUTH FOR DAMPING ID 



If damping is assumed unknown (but first mode frequency is known) , the distance 
measure curve becomes quite broad and well shaped. If the true system and the model 
agree in order, this measure reduces to zero at the true value. However, the lack 
of significance of distance measure indicates that the convergence of the MLE method 
will be poor and susceptible to noise. 
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.010 



.005 
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.000 



.002 .004 .006 

EstlMte of DMpIng Ratio 

Figure 12 
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.010 
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DISTANCE FROM 9 MODE TRUTH FOR DAMPING ID 

This conclusion remains if the true system is assumed larger than the model. In 
addition, a bias can be seen between the minimum distance damping ratio (the value 
the MLE method will converge to) and the true value. This bias is caused by the 
order mismatch and it is important to analyze this bias as it limits MLE performance. 
In addition, tradeoffs can now be performed to determine how large the identification 
model is required to be to keep the bias "small enough". 



49.260 



40.255 



40.250- 



40.245 



40.240 



40.235 




Estimate of Dapping Ratio 

Figure 13 
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CONCLUSIONS TO DATE 

We now have a tool which can be used to analyze these biases and other effects of 
system mismatch. 

The next step is to enlarge the study to include a more realistic situation. 

• Tool for Assessing Mismatch 

• Mismatch can result in biased estimate 

• Frequency estimate will need good initial guess 

• Damping Hard to estimate 

• More study needed 

Figure 14 

OUR PROGRAM 

We will be doing this by starting with a realistic strawman — the space shuttle 
with a payload on the end of the remote manipulator system (RMS) arm. In addition, 
we will be supporting and guiding the beam experiment as needed and will be guided 
by those results. 

Our first step is to determine representative parameters for sensor noises, actuator 
uncertainties, and other disturbances. This will help insure that our results are 
of more than academic interest. We will then be performing extensive identifiability 
studies to determine how accurately the parameters of the model can be identified. 

Finally, a controller will be designed which is compatible with the expected 
residual uncertainty and the performance compared to that obtained using a robust 
controller. 

t Use the Shuttle/RHS/Payload as stramman. 

- also beam experiment 
§ Determine Representative Environmental Parameters. 
• Investigate Identifiability of parameters. 
t Design a controller using available parameters. 
t Compare results with those using I LAS/Rate feedback. 

Figure 15 
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EXPECTED RESULTS 

In this way, we will determine the limitiations to be expected in performing on- 
orbit parameter identification. In addition, we will have investigated how these 
limitations affect the control design process and identified the achievable benefits 
using a high performance controller in place of a robust one. 

• Limitations on system identification 

• Effects of these limitations on control 

DESIGN. 

* how to design with available information 

* Benefits over 'Robust Control" 

Figure 16 
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EXPECTED NUMBER OF COMPONENT FAILURES IN A YEAR 

Why do we need to consider component unreliability? A large, lightweight structure in 
space will display many vibratory modes which may have to be actively damped to assure 
mission success. Effective control of these many modes will require use of a large 
number of sensors and actuators - possibly hundreds of them. Even if these control 
systems are serviced in orbit, one would like the service interval to be long - at 
least 1 year. With component mean time between failures which can reasonably be 
anticipated, one must expect many of the control system components to fail in the 
course of a year. Figure 1 presents the expected number of component failures per 
year versus component mean time between failures. 
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PROBABILITY OF ADDITIONAL FAILURES WITHIN 24 HOURS 



OF A COMPONENT FAILURE 



Will it be necessary to accommodate more than one component failure at a time? 
Following a component failure which is detected and identified, the system will be 
reconfigured to function with the remaining components . This reconfiguration should 
be accomplished in a relatively short time - perhaps on the order of 1 day. With the 
component mean time between failures which must be achieved to restrict Lhe expected 
number. of failures. in. a year to a reasonable value, the probability .of .one or^more 
additional failures in 1 day following a single failure is small. (See fig. 2.) To 
reduce this probability further, it may ;be possible, to '.reconfigure the failure 
detection and identification system more quickly than the complete control system. 




10,000 



100,000 
Component mean time between failures, br 



1,000,000 



Figure 2 



75 



SYSTEM OPERATIONAL APPROACH 
Figure 3 outlines the system operational approach. 



o 



Implement a quick-reaction standby mode. of control following 
indication of a component failure. 

Complete FDI system reconfiguration quickly— perhaps within 
a few hours, 

Complete control system reconfiguration as soon as possible— 
perhaps within a day. 

Figure 3 



SYSTEM DESIGN APPROACH 
Figure 4 outlines the system design approach. 



Utilize fault-tolerant computers as data processors for the controller 
AND FDI functions. 

Incorporate diagnostic and/or BITE information to give direct indication 
of component failures where possible. 

Implement an FDI system to detect and identify single additional control 

SYSTEM COMPONENT FAILURES— SENSORS AND ACTUATORS, 

Analyze specific system implementations to check for complete FDI coverage 

INCLUDING OTHER COMPONENTS SUCH AS POWER SUPPLIES, DATA CONVERTERS, DATA 
TRANSMISSION LINES, ETC, 

Figure 4 
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FAILURE DETECTION AND IDENTIFICATION 
Figure 5. outlines failure detection and identification. 



Many approaches to FDI have been suggested. 

Will not consider methods which require specification of the 
mode of failure, 

Will concentrate on methods which are applicable to both sensor 
and actuator failures. 

Two candidate techniques for initial consideration: 

An open loop method: Generalized parity relations 

A CLOSED LOOP METHOD: FAILURE DETECTION FILTER 

Figure 5 
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■'•■-.'.' GENERALIZED PARITY RELATIONS 

The usual concept of parity relations in FDI applies to redundant measurements such 
that a linear combination of the measurements can be set to zero. 

The generalized concept of parity relations utilizes measurements at more than one 
time, which introduces actuator input signals as well as sensor output signals. 



Discrete system model: 



x k+ i = A x k + B u k + n k 



£ k = C 2^ + D u k + v k 



Z k - C x k + D u k + v k 



Z k+1 -■eA^+^B.u ]c + ; D u k+1 + C r^ 



+, C n, ■ + v 



k+i 



H K+2 - CA * k + CAB H k + CB H k+1 + D H k+2 



+ CA n k + C n k+1 : 4- v^. 2 " 



etc. 














Example: 


















1 


0.3 


0.1 









A = 





' 1 - : 


0.4 


; B ■-■ 


0.1 









0.2 


0.5 




0.6 



c = 



10 



10 
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Except for the noise terms, 



,1k ~ A lk 



= x. 



2k A 2k 



*lk+l = X lk + °' 3 X 2k + °- 1 X 3k 

*2k+l = X 2k + °' 4 X 3k + °' 1 U k 

^lk+2 = X lk + °' 62 X 2k + °' 27 X 3k + °' 09 U k 

*2k+2 = 1 - 08 x 2k + °* 60 X 3k + °' 34 U k + °' 10 U k + 1 
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One simple generalized parity relation for a group of four rate gyros at one location 
can be derived considering: 

Measurements: y ik = c. T 0) k i = 1,2,3,4 ; 

The set of four cs in three-dimensional space are dependent, and one can find a. , 
such that . ■; ■.. 1 .■ 

4 
£ a. c. = 

.■ i=l .h 

Thus .,.-...■■ ■.'•: •• '...,' ... 

4 

r k = 2 a i y ik 

equals zero except for noise. A significantly nonzero r k indicates a failure among 
the four gyros ..: ... 

More generalized parity relations are formed from dependent rows of C, CA, CA 2 , 
etc. This eliminates the unknown state x k from the parity relations. 

An appropriate set of parity relations must be used, each depending on different sub- 
sets of sensors and actuators, so that from the set of residuals one can identify the 
failed component. ; 

One possible set of parity relations .is for the .example : 



r lk = 1 ' 9 ^Ik - °- 01 y 2k > 0.? Yl k+ i+ 0-9 y 2k+1 - y lk+2 
r 2k = -?.' 42 y 2 k +.. i: ? y 2 k + l - '.Yak+2 + °- 19 U k; + °' 1 u k + l 

r 3k = 4 yik + o.2 y 2k -- 4 y lk+1 * y 2k+1 - 0.1 u k 

Sensor 1 failure 
Sensor 2 failure 
Actuator failure 
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FAILURE DETECTION FILTER 

The failure detection filter (fig. 6) incorporates a model of the system being 
monitored. The difference between the filter output and the measured outputs of the 
system are fed back through the matrix D to cause the model to track the system. 
If the system is performing nominally, the output errors, after initial transients 
have subsided, are due only to unmodelled noises and disturbances and to model 
mismatch. 

When a component fails, the model no longer matches the behavior of the system and a 
significant output error develops. The filter gain matrix D is designed so that 
the output error due to each monitored failure is restricted to a single direction, 
and that direction is different for each component the filter is designed to 
monitor. Thus the presence of a significant output error - larger than that due 
to noise - indicates a component failure, and the direction of that error identifies 
the component that failed. 

It is possible to use detection filters to monitor sensor failures, actuator fail- 
ures, and to indicate significant mismatches in system parameters. 



Reference: Beard, R. V. 



Failure Accommodation in Linear Systems Through Self- 
Reorganization. NASA CR-118314, 1971. 
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QUICK-REACTION STANDBY CONTROL MODE FOLLOWING A FAILURE 
Figure 7 presents the quick-reaction standby control mode following a failure. 

a simple control system which maintains stability and restricts 
response to disturbances in spite of loss of an arbitrary component. 

Normal performance requirements may be suspended during standby operation. 

One candidate is a control system which only feeds back rate 
information to colocated actuators. 

Figure 7 
SYSTEM RECONFIGURATION 
An outline of the system reconfiguration is given in figure 8. 

Use of the interim standby mode allows some time in which to 
reconfigure the control system; a day seems reasonable. 

With the controllers implemented as digital data processors, 
reconfiguration amounts to revising the programs which these 
processors execute. 

This could be as difficult as designing the system in the 
first place — but using fewer sensors and/or actuators, 

it is assumed that most mission situations will permit 
ground-based support of the reconfiguration function with 
the results up-linked to the space system. 

Look for ways to simplify the reconfiguration process — perhaps 
using information generated during the original des-ign of 
the system. 

Figure 8 
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RESEARCH TASK 1 

Research Task 1 calls for development of a preliminary capability for the detection 
and identification of failures in the sensors and actuators of the laboratory 
experiment at Langley Research Center. 

la. Obtain from LaRC a dynamic model characterizing the laboratory apparatus 
and design a simple control system of some standard form to control the 
experimental plant. 

lb. Choose one or more of the currently available FDI techniques to apply to 
this problem. 

lc. Design the FDI processor for this application, including determination of 
all parameter values such as threshold settings, and check it out through 
simulation. 

Id. Prepare and deliver a FORTRAN program which can be compiled and executed on 
the 175 CYBER system at LaRC to demonstrate real-time FDI on the laboratory 
experiment system. 



RESEARCH TASK 2 

Research Task 2 investigates the possible advantages which may accrue from integrating 
into one information processor the functions of state estimation, sensor and actuator 
FDI, and parameter mismatch indication. 

2a. In the context of space structure control, separate the sensor and actuator 
FDI problem into two parts: one processor which monitors at least some of 
the sensors without use of system modeling, and one processor which monitors 
actuators and perhaps some sensors by a technique which depends on a model 
of the system dynamics. 

2b. Determine whether it would be efficient and convenient to extend the second 
FDI processor defined in Task 2a to monitor significant deviations in plant 
parameters from the modeled values. 

2c. Determine whether an additional price in computational capacity must be paid 
to derive an estimate of the system state from this . combined processor. 

2d. If the results of the above tasks have been favorable to the concept of the 
integrated processor, summarize the form that such a processor would likely 
take in a space structure control application and specify the logic which 
would be used to derive the best estimate of state under nominal and failure 
conditions . 
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RESEARCH TASK 3 

Research Task 3 studies the problem of efficient reconfiguration of a space structure 
control system, including its FDI processor, following a component failure. 

3a. Inquire into simple methods of reconfiguring the state estimator (s) follow- 
ing a sensor failure. Consider both the case of systems with highly 
redundant sets of sensors and the case of more sparse sensor sets. 

3b. Inquire into simple methods of reconfiguring the controller (s) following an 
actuator failure in the case of systems utilizing only linear control 
devices. Consider both the case of systems with highly redundant sets of 
actuators and the case of more sparse actuator sets. 

3c. Inquire into simple methods of reconfiguring the controller (s) following an 
actuator failure in the case of systems utilizing on-off jets. Consider 
both the case of systems with highly redundant sets of jets and the case of 
a more sparse complement of jets. 

3d. Inquire into simple methods of reconfiguring the FDI processor following a 
sensor or actuator failure. 



RESEARCH TASK 4 

Research Task 4 investigates design techniques for multidimensional systems utilizing 
on-off control devices. 

4a. Determine the nature of optimal trajectories for the problem stated in the 
Technical discussion in the case of a single-input system. 

4b. Define the nature of the cost function J(x) for the single-input system for 
initial conditions near the origin and far from the origin. Try to charac- 
terize this function on the basis of analytic arguments and check the results 
with computed optimal trajectories. 

4c. Synthesize the controller which minimizes the Hamiltonian function with the 
costate variable identified with the gradient of an approximate cost func- 
tion. Compare the performance of the resulting feedback control system 
with optimal performance . 

4d. Attempt to generalize the above steps to multidimensional systems with an 
arbitrary number of on-off actuators. 
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INTRODUCTION 

The Jet Propulsion Laboratory flexible beam facility was designed to be sufficiently 
general to demonstrate a variety of flexible body control concepts. Some of these 
concepts include static shape control, vibration control, multi-input/multi-output 
control, noncolocated sensors and actuators, adaptive control , and distributed con- 
trol. The facility consists of a 3.8-m by 15.2-cm by 0.08-cm (12^-ft by 6-in by 
1/32-in) pinned-free flexible stainless-steel beam hanging in the Earth's gravity. 
Four eddy- current position sensors and three brushless dc motors provide position 
information and force capability. A 6502 microprocessor with developed software 
completes the control loop through twelve bit D/A and A/D interfaces. 

• DEMONSTRATE LSS CONTROL IDEAS 

• DEVELOP FACILITY WITH GENERAL APPLICABILITY 



STUDY HARDWARE LIMITATIONS 



HINGED -FREE 
FLEXIBLE BEAM 




JPL 
CELESTARIUM 



CRANE CARRIAGE 
FOR HOISTING/ 
LOWERING OF 
BEAM/TOWER 



.6.1 m 
(20') 



Figure 1 
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BEAM DYNAMICS - STRING OR ELASTIC BEAM? 



The hanging beam is in tension due to. gravity, and because of the very small bending 
stiffness, low- frequency modes appear string- like, while the higher frequency modes 
appear beam-like. Dynamic mode shapes are computed using string or beam partial 
differential equations of motion, and the results are then combined to produce accu- 
rate system mode frequencies. 
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MODES 



n =1 
0. 30 Hz 
(RIGID 
BODY) 




n = 2 
0. 74 Hz 





n =4 
2. 00 Hz 





F igujre 3 
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SYSTEM NONLINEAR! TIES 



The effects of all major system nonlinearities are analyzed using an open-loop 
perturbation analysis to determine the frequency shift or damping ratio shift due 
to the perturbing force. The largest effect is due to aerodynamic drag with the 
result being the addition of 1% of critical damping. The remaining list of non- 
linearities were constrained to have a negligible impact on the beam dynamics 
from the beginning of the design process using the perturbation analysis to pre- 
dict the resulting magnitudes of the various nonlinearities . 



• AERODYNAMIC DRAG € q q 

• FINITE TRANSDUCER SIZE 



/ y(x) dx 



• ANGLE DEPENDENT TORQUE T » T c (1 - cq) 

• NON LINEAR TORQUE T = e6 3 

• NON LINEAR DAMPING T = e6 3 

• ADDITIONAL MASS/STIFFNESS p(x)/EI(x) 

OPEN-LOOP PERTURBATION ANALYSIS 

p 
q + cj q = € f (q, q) 

Figure 4 
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CONTROL-LAW DESIGN 



For the initial studies into the problem of vibration control, optimal control theory 
was used to produce estimators and control laws. All results were obtained using 
discrete time design techniques. 



LINEAR QUADRATIC ESTIMATOR AND REGULATOR 
x(k + ll = $ x(k) + G u(k) +rw(k) 
z(k) = H x(k) + v(k) 

CONTROL MIN X) * Ax + u T Bu s.t. x(k + l) = $x(k) + Gu(k) 

ESTIMATOR MInJ]w T Q" 1 w + v T R _1 v s.t. x(k + l) - «x(k) + Gu(k) 
RESULTS (DOPTSYS) u = C x 

x~—&x + Gu + K (z - Hx) 

Figure 5 
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MICROPROCESSOR IMPLEMENTATION 



Software for a 6502 microprocessor was developed to implement any linear, time 
invariant estimator and controller, and was placed into IK of ROM. IK of RAM 
is used to store the particular control system matrices and the scratch space 
for the ROM software. The software samples the sensor vector, updates the state- 
vector estimate, and outputs the control vector. The loop time is dependent on 
the dimensions of these vectors with a six-state controller using 1 sensor and 
1 actuator requiring 50 msec per update. Twelve-bit resolution is retained in 
the D/A and A/D processes, while sixteen bits are used for internal computation. 
Fixed-point double-precision arithmetic is performed entirely in software. 
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RESULTS 



These preliminary results show the free end response of the beam to various inputs, 
The open-loop initial condition (IC) and impulse responses are shown. Next, low- 
gain closed- loop responses to both of the above inputs are shown. The controller 
in this case consists of a single-position sensor and a single-force actuator at 
the free end. Using a six-state controller, about 40 percent of critical damping 
is achieved. The small residual is at the frequency of the first unmodeled mode. 
By decreasing the cost of control in order to obtain improved performance, the 
first unmodeled mode is destabilized. 

Shape control, distributed control, and refining of the parameters describing the 
beam facility are all activities planned for the future. 
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BASIC ELEMENTS 

Parameter Identification Subsystem - 

Provides "local" estimate of system model parameters at a specific operating condi- 
tion. 

Techniques that have been used : 

Newton-Raphson parameter identification 
Output-error parameter identification 

Parameter Extrapolation (Learning) and Memory Subsystem - 

Provides "local" extrapolation of system parameters as a function of measured con- 
figuration variables. 

Technique that has been used: 

Least-squares fit of functional model relating model parameters to configuration 
variables. 

Controller - 

Uses functional model stored in Memory Subsystem to generate control forces. 

PARAMETER ADAPTIVE CONTROL SYSTEM 
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APPLICATION OF ADAPTIVE/LEARNING 
APPROACH TO DISTRIBUTED SYSTEMS 
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FEATURES OF ADAPTIVE/LEARNING APPROACH 



Approach combines advantages of a priori modeling and gain scheduling with those of 
on-line identification and adaptive gain computation. 

Convergence criteria may be used to monitor performance of parameter-identification 
and parameter-extrapolation subsystems. 
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APPLICATION OF ADAPTIVE/LEARNING CONTROL TO 
DISTRIBUTED SYSTEMS - PARAMETER IDENTIFICATION 

Response Representation - 

NM 
y(x,t,£) =2 ?. (x,A)w. (t,p (£)) 
i=l 1 1 

£ . - spatial approximation function 

w. - modal amplitude approximation 

p. - parameter vector in difference equation model for the evolution of w 

I - configuration variable 

Model Parameter Identification - 

Adjust p. to obtain best fit of model data to sensor data. 

w. (k) = A, .w. (k - 1) + A_.w. (k - 2) + B F (k - 1) + B F (k - 2) 
1 ll 1 Zl 1 11 1 A2. X 

Pi = [ A li' A 2i' B li' B 2i] 

F. - force command input 
i 

Approximation Function Tuning - 

NB 

h*' l) =2 c ik i k (x '« 

k=l 

a (x,£) - a set of basis functions 

ik 
C - response model spatial function tuning coefficients 

ik 

STRUCTURE OF CURRENT IDENTIFICATION PROGRAM 

Actual mode shape information 
for calculating H . 
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w(k) - vector o£ modal amplitude estimates 

p - vector of model parameters 

V (k) - vector of displacement measurements 
m 



v(k) - error vector 

H . - matrix related to assumed 
- 1 spatial functions ? 
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SIMULATION STUDIES OF FLEXIBLE BEAM UNFORCED RESPONSE - SPATIAL VARIABLE PLOTS 
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PROPOSED BEAM-EXPERIMENT EXTENSION TO TEST - VARIABLE-CONFIGURATION STRUCTURE 




Variable-configuration experiment 




Variable attachment point- 



Beam 



1 



Configuration at test 1 



U JL 



U jIl 2 



Configuration at test 2 



Figure 7 
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PROBLEMS FOR INVESTIGATION 
L. Development of function tuning algorithm. 

2. Development of test-initiation algorithm. 

3. Evaluation of state-estimation and parameter-identification algorithms. 

4. Development of variable-configuration experiment. 

5. Performance evaluation of closed-loop adaptive/learning control system. 
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OBJECTIVES FOR EXPERIMENTAL FACILITY 

1. Very low structural damping, ? < 0.001 

2. Noncolocated sensors/actuators 

3. Changing (and uncertain) inertias 

SKETCH OF TORQUER/ENCODER MECHANISM 



/// // , . 
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TRANSFER FUNCTIONS 
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NONCOLOCATED CONTROL 
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NONCOLOCATED CONTROL WITH APPENDAGE 
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STATEMENT OF TASK 

At NASA's request, the Aeronautics and Space Engineering Board established an 
ad hoc committee to address the following tasks: 

1. Evaluate and provide comments and recomendations on the objectives, ap- 
proach, content, and technical balance of NASA's Large Space Systems Tech- 
nology Program plan. 

2. Review and recommend means to assure an effective exchange of information 
on technological developments in the area of large space systems between 
and among NASA, its contractors, and user industries and agencies. 
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TRIANGULAR STRONGBACK-MEDIUM PLATFORM 
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GENERAL CONCLUSIONS AND RECOMENDATIONS 

1. LSS missions furnish solid justification for continuation of LSST program. 
"Current Plan" inadequate for fulfilment of needs during next 5 years. 

2. LSS program should be integrated, with designated Headquarters offices as over- 
all director. 

3. Prioritization of mission model essential to focus the technology program. 

4. Flight experiments during next 5 years essential to focus TLSS program. 

5. "Roadmaps" needed for technological development: 

a) Total system and 

b) Local for each technology discipline. 

6. Study sensitivity of system performance to relaxed requirements on technology. 

7. Large platform with several experiments a particularly imaginative concept; 
further study needed. Can products of different users be combined ? 

8. Man at GEO needs to be vigorously examined. Committee not convinced benefits 
will outweigh costs. 

9. Transfer of technology and coordination of work with other agencies (notably 
DOD) quite adequate and improving. 
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TECHNOLOGICAL CONCLUSIONS AND RECOMMENDATIONS 

1. Reinforce certain areas in TLSS 5-year plan: 

a) Figure determination and control 

b) EM interference 

c) Automated vs. manual assembly 

d) Docking 

e) Rel., repair, and maintain 

2. Comparative evaluation of alternatives required. 

3. Definition of orbital transfer vehicles lacking. Electric propulsion should be 
pursued, but capability in chemical should be preserved. 

4. LSST plan tends to emphasize structures more than certain other vital areas. 

5. JPL study of ODSRS should be evaluated; considered a rather unlikely candidate. 

6. SEP has a significant role for orbital transfer, station-keeping, and planetary 
missions. '-'■-.' 

7. More attention needed to locations and systems for data processing. Needs a 
significant start. 

8. Traffic models must be better developed to define future COMSATS. 

9. Efforts applauded on low-cost ground stations. Success key is a large market. 

10. Current program reasonable relative to large-platform COMSATS at GEO. 

11. In structures and materials, emphasize optimization, loads determination, dy- 
namics, damping, integrated design, and special materials for space. 

12. Activity on deployment, assembly, and fabrication is effective. Time for eval- 
uation of alternative approaches . 

13. Technology of active figure sensing and control needs greater attention. 

14. LSS human factors program endorsed, but specific line item should be added 
within LSST budget. 
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OBJECTIVES 

Generic model of a planar space structure having infinite number degrees of freedom 
(DOF) : 3 identical, simple beams joined with ball-socket connectors. 

Demonstrate roll, pitch, and yaw attitude stabilization using angle sensors and tor- 
quers at the midpoints of all 3 beams . 

Find exact transcendental transfer functions from torques to angles producing exact 
poles and zeros . 

Many identical natural frequencies in roll and pitch cause no control problems. 

Compare closed-loop predictions of truncated models to infinite DOF model. 

Compare finite-element approximate models with exact infinite-element model. 
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NOMENCLATURE AND UNITS 



a = mass/unit length of beams 
EI = bending stiffness of beams 
21 = length of beams 
y = out-of-plane displacement 
6 = out-of-plane slope 
Q = in-place bending moment 
F = out-of-plane shear force 
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TRANSFER MATRICES 
From partial differential equation of simple, we can obtain T(X) such that: 



8 eqs. 



where 



T(X) = 



x (X) = T(X)x (X) 
x 4 (X) = T(2X)x 3 (X) 



T T T T 
1 2 3 4 

T T T T 
4 1 2 3 

T T T T 
3 4 1 l 2 

T T T T 
L 2 3 4 1 



= —(cosh A ± cos A) 



= — (sinh A + sin A) 



For roll motion: 
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(antisymmetry) 
(ball-socket joints) 
(compatibility of displacements) 
(transmission of shear force) 
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ROLL TRANSFER FUNCTION 
Combining 8 boundary conditions with 8 equations of two transfer matrices gives 



<f>(s)/Q R (s) 



where 



Q R (s) 



(s) = EIX N(s) - 



D(s) 



(dimensional) 



N(s) = Mcosh 2X sin 2X - sinh 2X cos 2X) (cosh X sin X - sinh X cos X) 

-h sinh 2X sin 2X(1 + cosh A cos X) 
D(s) - (cosh 2X sin ! 2X - sinh, 2X cos 2X) (sinh^.X sin X + sinh 2X sin 2X) 

+ sinh 2X sin 2X (cosh X sin X - sinh X cos X) 
We may now write: 



(s) 



Q R (8) 



2 . 2 
1 °° 1 + s/Z , 

~2 n " 2~2" 

Js. n=l 1 + s ./a) . 



(dimensional) 



where 



J = 2ai = Roll moment inertia of rigid truss 



Using a root solver code, we can find exact poles and zeros (w ,Z ) from tran- 
scendental transfer function above. First 6 deformation modes and rigid-body mode are 
shown in figure 2. Zeros are for colocated angle sensor. (Note that poles and zeros.' 
alternate ! ) 



o ZERO 
x POLE 



x 47.012 



6 4 1 .322 
* 39.478 



25.081 
24 . 249 



18.135 



x 13.799 

o 10.846 
x 9.869 



3.747 
3.162 
I .586 
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ROLL-MODE SHAPES 
From the partial differential equation of the simple beam, 

y(n) = y^ttn) + e^un) - Q^tXn) + F^Un) (Beam i, o<n<i) 

y(C) = -y^US) - e 4 T 2 (A?) + Q 4 T 3 (A5) - F^U?) (Beam II, < £ < 1) 
For antisymmetric roll motion, Y;L = Y/ ^ = ^ = Q ^ = . For each co yielding X , 
find 6 lf 9 4 , F 1 , F 4 from transfer matrices ("eigenvector") derived from mode 
shapes from above. The first 6 mode shapes are shown in figure 3. 

ROLL-RATE STABILIZATION 

Colocate roll-rate sensor with roll torquer. Use simple proportional feedback 

Q R = -K<j>. Figure 4 is part of the infinite-dimensional root locus vs. K. All modes 

are stabilized: Only 2 deformation modes are significant. Figure 5 is the root locus 
vs. K for the truncated model including only the 2 significant deformation modes. 
Figure 6 shows the 6th order model in more detail. 
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ROLL-ATTITUDE STABILIZATION 



Colocate roll-angle sensor with roll torquer. Use lead compensation 

s + 1 



-K 



'(s) 



*R s + 10 
Figure 7 is the root locus vs. K for the infinite-dimensional model. 



PITCH CONTROL 

Figure 8 shows the free-body diagram for analyzing pitch motions using 2 torquers 
Qp and Q p • Using the same method, we found the exact transfer function 6(s)/Q (s) 

1 2 p 

The first few poles and zeros are shown in figure 9 for colocated sensors and torquers. 
Note the pole/zero cancellation plunge modes. The first few mode shapes are shown in 
figure 10. Some, but not all, frequencies are identical to the roll frequencies. 
Sensors do not pick up roll motions so there is no roll-pitch coupling. 
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YAW CONTROL 

Figure 11 shows the free-body diagram for analyzing yaw motions by using 3 tor- 
quers (at beam midpoints) . We must now include the axial deformation of the beams 
since bending is in-plane. The values (x.,G.) are the axial displacements and axial 
forces at the points indicated. 

Using the same methods , we found the first few exact poles and zeros for co- 
located sensors/torquers , shown in figure 12 for 3 different beam-slenderness ratios. 
The first few mode shapes are shown in figure 13 . 
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CONCLUSIONS 

1. Exact infinite degree-of -freedom transfer functions can be found for truss/frame 
space structures using simple beam "transfer matrices". 

2. By using symmetries, we can omit unexcitable modes in numerator/denominator, 
saving computation time. 

3. Proper location of sensors and actuators avoids feedback-coupling of naturally un- 
coupled modes at the same resonant frequency. 

4. Transfer function zeros (which depend on mode shapes) are found more accurately by 
this method than by the usual "finite-element" codes. 
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CONSIDERATIONS FOR DYNAMIC SIMULATION 
OF LARGE SPACE STRUCTURES 

- Large attitude changes 

- Small deformations 

- Arbitrary initial conditions 

STRUCTURAL DYNAMIC ISSUES IN LARGE SPACE STRUCTURES 

Lump masses at nodes ? 

Treat members as linear springs ? 

Formulate the differential equations governing the coordinates of the nodal parti- 
cles ? 

Solve the differential equations ? 

N nodes (N very large) 
3N second-order differential equation 
6N third-order differential equation 

The fact that only small deformations are of interest is not exploited. 

Small motions can be described in terms of principal modes - and only a few modes 
suffice. Can one use the same approach to describe large motions ? (See fig. 1.) 





Figure 1 
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SIMULATION CONFIGURATION 



X x X i 5 >\ 



*i>*+ y 




Configuration Variables 
F x. 



V 



or 



Y l Y 6 ; V Z 2 

Modal Coordinates 



Figure 2 



*v 



v 



y . = x. + £ A q 
D J k=l : k k 

Generalized Speeds 



u„ = 



' U 2 = ^%' U 3 = ^-%>' U 3+k = \ (k 



= 1,..., v) 
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EQUATIONS OF MOTION 



Kinematical 



= u. 



z = u cos - u sin 
z = u sin 6 + u cos 6 

\ = U 3 + k (k = 1 V) 

Dynamical 

v 3 

Xi l " w 2*2 + W 3^3 + 2 2 m i (Y 2i-l A 2ij ~ ^Ai-j^+j 

j=l i=l 

3 r ' l 

= 2 t Q 2i " m i a i2 )Y 2i-l " (Q 2i-l - m i a il )Y 2l J " Vl + ^ 
1=1 

3 3 

= 2 <e 2i _ 1 - v±i> vi + ™ 3 = 2 (Q 2i - m i a i2 ) 

x=l i=3 



3 3 

2 m i (y 2i A 2i-k:.- y 2i-l A 2ik )G l " *3 + k = - G k + \\ + X 
1=1 i=l 



^ m (v A -v A )u - u =-G+X, q+5^m. (A„, , a. + A_ ., a. „) 
ft* i Y 2i 2i-k: Y 2i-1 2ik 1 3+k k k H k f* i 2l ~ k i ^ik i2 

(k = 1 , . . . , v) 



Initial Conditions 

Given: x (0) , ^(0), x ? (0); ^(0), x (0) , x 3 <0) 

Needed: z (0) , % 2 (0) , 0(0), ^(0), ^(0), ^(0); q fc (0), u 3+k (°) < k = !/■••/ v) 

6 + 2v unknowns 
Equations : 

Zl (0) - x 2 . (0)0(0) + V A 2i _ k q k (0) = x 2 ._ r (0) - x 2 ._ ± 



k=l 

z 2 (0) + x 2± (0)9(0) + f, A 2ik q k (0) - x 2 . (0) - x 2 . 

k=l 

-y 2± (0)u (0) + u (0) + 2 A 2 ._ k u 3+k (0) = x 2 ._ l( 0) - x 2 . (0)0(0) 

k=l 

- y 2i-l ( ° )U l (0) + U 3 (0) + £ A 2ik U 3 + k (0) = - k 2i-l { ° )Q{0) + *2i (0) 

k=l 

(i = 1, 2, 3) 

12 equations 
If v < 3, the number of equations is greater than the number of unknowns. 
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SIMULATION RESULTS 



5~0b0 A/ 




MOOO N 



3000 A/ 



Figure 3 



US', 60 




Figure 4 (<j) = exact; $ = linear solution) 
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